Majority-vote model on random graphs 
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The majority- vote model with noise on Erdos-Renyi's random graphs has been studied. Monte 
Carlo simulations were performed to characterize the order-disorder phase transition appearing in 
the system. We found that the value of the critical noise parameter qc is an increasing function of the 
mean connectivity z of the random graph. The critical exponents P/v, ■y/v and l/v were calculated 
for several values of z, and our analysis yielded critical exponents satisfying the hyperscaling relation 
with effective dimensionality equal to unity. 
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I. INTRODUCTION 

Since Erdos and Renyi's work more than forty years 
ago pi, ^] , intense theoretical research on random graphs 
has been taking place 0, 3 . In particular, models of net- 
works with more complex connectivities than the tradi- 
tional uncorrelated random graphs have been introduced 

1^ and used to describe many systems in nature and 
society 011 0. 

A random graph is a set of N vertices (sites) connected 
by B links (bonds). The probability p that a given pair 
of sites is connected by a bond is p = 2B/N{N — 1). 
The connectivity of a site is defined as the total num- 
ber of bonds connected to it, that is ki = Cij, where 
Cij = 1 if there is a link between the sites i and j and 
Cij = otherwise. Random graphs are completely char- 
acterized by the mean number of bonds per site, or the 
average connectivity z — p{N — 1). In the limit N ^ oo 
the distribution of connectivities is given by the Poisson 
distribution. 

For values of z < 1 the random graph does not have a 
percolating cluster 0, ^ . There are a few disconnected 
clusters and there is no long-range order on such systems. 
For 1 < z < 4 there is a percolating cluster, but there 
are a few small islands disconnected from the giant clus- 
ter. These small islands do not contribute to the system 
dynamics and so they are excluded from our simulations. 
For values of z > 4 almost all the sites belong to the 
giant cluster, so no sites need to be excluded from the 
dynamics. 

Our goal in this work is to identify the critical char- 
acter of the majority-vote model with noise on random 
graphs. Previous works on the majority-vote model con- 
sider the spins either on regular d-dimensional lattices 
[TTl IT^ or on small- world networks 13, 14] interpolat- 
ing between regular lattices and random graphs. We use 
Monte Carlo (MC) simulations and standard finite-size 



scaling techniques to determine the critical noise param- 
eter qc, as well as the exponents P/v, ^/v and l/v for 
several values of the mean connectivity z of the graph. 
We also use mean-field approximation to obtain the phase 
diagram in the qc — z space and make a comparison with 
the corresponding phase diagram that follows from our 
simulations. 

This paper is organized in the following way: In sec- 
tion^we describe the isotropic majority- vote model with 
noise and introduce the relevant quantities used in our 
simulations. In section Ulll we present our results along 
with a discussion. And finally, in the last section we 
present our conclusions. 



II. THE MODEL AND COMPUTATIONAL 
METHOD 



The isotropic majority-vote model with noise is defined 
by a set of spin variables {ui}, where each spin is asso- 
ciated to one vertex of the random graph and can have 
the values ±1. The system dynamics is as follows: For 
each spin we determine the sign of the majority of its 
neighboring spins, that is, all the spins that are linked to 
it. With probability q the spin takes the opposite sign of 
the majority of its neighbors, and it takes the same sign 
with probability (1 — q). The probability q is known as 
the noise parameter. 

The probability of a single-spin-flip is given by 



w{(yi) = 2 



l-(l-2(?)a,5(^a,+5) 



5=1 



(1) 
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where S{x) — sgn(a;) if x 7^ and 5(0) = 0, and the 
summation is over all the ki spins connected to the spin 
at site i. 

To study the critical behavior of the model we consider 
the magnetization Mm, the susceptibility xn, and the 
Binder's fourth-order cumulant Un . These quantities are 
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defined by 



N 



Tl C 



XN{q)=N[{ {m')T)c - { {m)T)h] 



UNiq) - 1 
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where N is tlie number of vertices of the random graph 
with fixed z, (...)t denotes time averages taken in the 
stationary regime, and (...)c stands for configurational 
averages. 

These quantities are functions of the noise parameter q 
and, in the critical region, satisfy the following finite-size 
scaling relations |3 



MN{q) ^N-P'''M{N^'''e) 



(5) 
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FIG. 1: Magnetization and susceptibility as a function of the 
noise parameter g, for TV = 10000 sites. From left to right we 
have 2 = 2, 4, 6, 8, 10, 20 and 50 



III. RESULTS AND DISCUSSION 



XN{q) = N->^''x{N'/-e) (6) 

UN{q) = U{N'/''e) (7) 

where e = q — qc- From the size dependence of Mjv and 
Xn we obtained the exponents /S/v and 7/;^, respectively. 
The maximum value of the susceptibility also scales as 
N'y/'^ _ Moreover, the value of q for which xn has a max- 
imum, qc{N), is expected to scale with the system size 
as 

g,(iV) = g, + 5iV-i/^ (8) 

where the constant b is close to unity. The above relation 
is used to determine the exponent l/iy and also to provide 
a check for the values of qc obtained from the analysis of 
the Binder's cumulant (Eq. [TJ. Finally, we have checked 
whether the calculated exponents satisfy the hyperscaling 
hypothesis 

2l3/u + ^lu = D,ff (9) 

in order to get the effective dimensionality, D^f / , for sev- 
eral values of z. 

We have performed Monte Carlo simulations on ran- 
dom graphs with various values of mean connectiv- 
ity z. For a given z, we used systems of size N = 
1000, 1750, 2500, 5000 and 10000. We waited 8000 Monte 
Carlo steps (MCS) to make the system reach the steady 
state, and the time averages were estimated from the next 
4000 MCS. In our simulations, one MCS is accomplished 
after all the N spins are updated. The simulations were 
performed using the standard C random number gener- 
ator. For all sets of parameters, we have generated ten 
distinct random networks, and we have simulated ten in- 
dependent runs for each distinct graph. 



In Fig. n we show the dependence of the magnetiza- 
tion and the susceptibility xn on the noise param- 
eter, obtained from simulations on random graphs with 
N = 10000 nodes and several values of the average con- 
nectivity z. In part (a) each curve for M^r, for a given 
value of N and z, clearly indicates that there exists a 
phase transition from an ordered state to a disordered 
state. We also notice that the transition occurs at a value 
of the critical noise parameter gc, which is an increasing 
function of the mean connectivity z of the random graph. 
In part (b) we show the corresponding behavior of suscep- 
tibilities Xn ■ The value of q where xn has a maximum is 
here identified as qc{N). In Fig. |2we plot the Binder's 
fourth-order cumulant Un for different values of N and 
two distinct values of z. The critical noise paramenter q^ 
for a given value of z, is estimated as the point where the 
curves for different system sizes N intercept each other. 
In this way we have obtained the phase diagram shown 
in Fig. El 

The phase diagram of the majority-vote model on ran- 
dom graphs shows that for a given graph (fixed z) the 
system becomes ordered for q < q^ whereas it has zero 
magnetization for q > qc- We notice that the increase of 
qc is more pronounced for small values of z. The error 
bars in qc (see Tabled are smaller than the symbols. In 
the figure, it is also shown the values of qc obtained from 
mean-field (MF) approximation. For small connectivi- 
ties z the MF estimate of the critical noise parameter is 
very inaccurate. In particular MF theory gives qc = 
for z <2, whereas our MC phase diagram exhibits an or- 
dered state for all values of the mean connectivity greater 
than one. This is in agreement with the limiting value of 
2: = 1 for the existence of a percolating cluster and, there- 
fore, the onset of long-range order in the system. How- 
ever, as z increases the two estimates get closer. This is 
expected because MF approximation becomes more pre- 
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FIG. 2: Binder's fourth-order cumulant as a function of q 
and five values of system size A''. In part (a) we have z = 8 
and in part (b) z = 20. 
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FIG. 4: ln(Mjv(ijc)) versus InA'^. From top to bottom, z 
2,4,6,10,20,50,100. 
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FIG. 5: (a) Plot of ln(xjv(gc)) and ln(xjv(max)) versus In A'^. 
(b) The dependence of the noise parameter qc{N) on system 
size. The extrapolation gives an independent estimation for 
Qc- The data are for the case of mean connectivity z — 3. 




FIG. 3: The phase diagram, showing the dependence of the 
critical noise paramenter qc on the average connectivity z, 
obtained from MC simulations and from MF approximation. 



cise as the number of interacting nodes is increased. We 
have also performed simulations for random graphs with 
higher values of z, such as 2 = 50, z = 100, and z = 1000. 
The corresponding values of the critical noise (not shown 
in the figure) are smaller than 0.50, which is the limiting 
prediction value as provided by mean-field theory when 
z — > 00. 

In Fig. ^we plot the dependence of the magnetization 
at q = qc with the system size. The straight lines, ob- 
tained from simulations with different values of the mean 
connectivity z, confirm the scaling of the magnetization 
according to Eq. The slopes of curves correspond 

to the exponent ratio P/v. Our results show that the 



increase of P/v with z is quite small. 

In Fig. EJa) we display the scalings for the susceptibil- 
ity at q = qc, XNilc), a-iid for its maximum amplitude, 
XAr(max). The exponent ratio j /v are obtained from the 
slopes of the straigth lines. For almost all the values of 2, 
the exponents 7/1/ of the two estimates agree within error 
bars (Table P). An increased 2 means a slight tendency 
to decrease the exponent ratio j/i'. 

In a similar way, for fixed z the critical exponent l/v 
was obtained from a plot of In qc{N) — qc versus In N (see 
Eq. ©). We used the corresponding values of qc{N) 
that follow from the maximum of the susceptibility and 
the limiting value of qc which has been obtained from 
Binder's cumulant. The slope of the resulting straight 
line equals the exponent 1/v. The results quoted in Table 
U indicate that l/i^ is not a monotonic function of the 
mean connectivity z. 

Fig. [S{b) illustrates the scaling relation for qc{N) given 
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FIG. 6: Data-collapsing for five different values of A'', with 
z = 10. 



in Eq. JSJ. The constant 6 equals the slope of the straigth 
line, whereas the extrapolation of the fitting provides an 
alternative way of determining the critical parameter qc. 
We have obtained a quite satisfactory agreement between 
the values of qc determined in this way and the corre- 
sponding ones that follow from the analysis of Binder's 
cumulant. 

In Fig. El we show the data-collapse plot for M(u) ~ 
Miq{q)N^^'^ , which is a universal function of the com- 
bined variable u — N^/^{q — qc). We have also obtained 
quite good data-collapse for x{u) = XN{q)N~'^/^ . The 
collapsing of curves for five different system sizes corrob- 
orates the quoted values for qc, 13 /v, ^/v and 

Table m resumes the values (along with errors) of gc, 
the three critical exponents {-^/v was obtained by using 
two different scalings), and the effective dimensionality 
of the system. It is worth to mention that, for all z, 
the value Dcff = 1, which has been obtained from the 
hyperscaling hypothesis (Eq. 0), is satisfied when we 
use both estimate procedures for the exponent ratio ^/v. 

As far as we know, there is no previous works study- 
ing the majority-vote model on Erdos-Renyi's graphs, to 
allow a direct comparison of the present results. Yet, for 
completeness, it would be of interest to mention earlier 
simulations of the maiority-vote model on other kinds of 
networks. Campos et al [T^ investigated the phase di- 
agram and critical behavior of the majority-vote model 
on small-world networks |^ by rewiring the two dimen- 
sional square lattice. Using a similar procedure to ours 
they found critical exponents depending on the fraction 
of long-range interactions and satisfying the hiperscaling 
relation with Dcff = 2 (the dimensionality of the reg- 
ular lattice). On the other hand, the model which has 
been defined on a regular lattice has critical exponents 



that fall into the same class of universality as the corre- 
sponding equilibrium Ising model [Til hj. The results 
of the present simulations show that the majority- vote 
models defined on a regular lattice, on small-world net- 
works, and on Erdos-Renyi's random graphs belong to 

TABLE L The critical noise qc, the critical exponents, and 
the effective dimensionality Dejf, for random networks with 
mean connectivity z. 



P/v 'y/u 



l/u D. 



2 

3 

4 

6 

8 
10 
20 
50 
100 



0.066(1) 
0.1349(6) 
0.181(1) 
0.2403(5) 
0.2753(3) 
0.2998(4) 
0.3586(2) 
0.4110(2) 
0.4368(3) 



0.24(2) 
0.233(2) 
0.242(6) 
0.245(7) 
0.242(5) 
0.259(1) 
0.255(4) 
0.271(4) 
0.267(4) 



0.51(3) 
0.529(2) 
0.54(1) 
0.514(4) 
0.510(9) 
0.483(5) 
0.501(6) 
0.465(6) 
0.467(7) 



0.52(2) 
0.529(7) 
0.515(6) 
0.507(4) 
0.510(3) 
0.502(5) 
0.503(2) 
0.485(4) 
0.479(4) 



0.48(5 
0.37(9 
0.59(7 
0.44(6 
0.56(3 
0.51(3 
0.49(4 
0.39(5 
0.47(3 




"Obtained using Xiv(9c)- See Eq. l6t . 
''Obtained using xjvC""^^:). 
'^Estimated using j/u from Xn{Qc). 



different universality classes. 

Finally, we remark that our MC results are quite dif- 
ferent from the mean-field estimates j3/v = 1, ^/v = 2 
and \/v = 2, which result in Z? = 4 for the upper critical 
dimensionality. This is a reasonable result since for all 
networks simulated we are far away from the mean-field 
picture where every spin interacts with all the remaining 

— 1 spins. 



IV. CONCLUSION 

We have obtained the phase diagram and critical expo- 
nents of the majority-vote model with noise on random 
graphs. The second-order phase transition which occurs 
in the model with mean connectivity z > 1 has exponents 
that show a slight variation along the critical line. Never- 
theless, our Monte Carlo simulations have demonstrated 
that the effective dimensionality Dcf / equals unity, for all 
values of z. This interesting result may suggest that other 
spin models defined on random graphs have exponents 
which satisfy the hyperscaling relation with D^ff = 1. 
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